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_ Abstract 

o 

I This paper provides a full geometric development of a new technique called un-reduction, for dealing with dynam- 

ics and optimal control problems posed on spaces that are unwieldy for numerical implementation. The technique, 
^N) which was originally concieved for an application to image dynamics, uses Lagrangian reduction by symmetry in 

^ reverse. A deeper understanding of un-reduction leads to new developments in image matching which serve to 

^ illustrate the mathematical power of the technique. 

Q 

^ 1 Introduction 

Recently there has been a certain degree of interest in geodesic shape matching, as in [ , , ] and references 
therein. This interest has grown from the desire to assist applications in medical imaging by providing a coherent, 
^ quantitative method for comparing shapes. 

For our purposes, a shape is taken to be an embedded or immersed submanifold of an ambient space. For example, 
Ch a closed, simple, planar curve is a shape since it is the image of an embedding ^ M^. The space of shapes 

(or shape space) may be realised as a quotient. Given two smooth manifolds Ai and Af, the shape space A^-type 
T— I submanifolds of Af may be identified with the quotient space E := Emb {A4,Af) / Diff {A4). That is, the space of 
J> embeddings of A4 into Af up to re-parametrization of A4. 

In order to compare two shapes, a natural way to proceed is to construct a path between them. This is known 
as the matching problem. For example, the matching problem for closed, simple planar curves may be stated as 

O follows: Let po and pi be two smooth submanifolds in of 5^-type. Find a path of submanifolds of S'^-type, p{t), 
such that p{0) = po and p(l) = pi. In practice there are many such paths, therefore a variational principle is used 

T-H to select a preferred path. 

o 

^ Problem 1.1 (The matching problem) Minimize the functional 



X 



{p,p) dt 



subject to the boundary conditions p{0) = pa and p(l) = pi, where E denotes the space of all submanifolds of -type in the 
plane and P : TE — > M fs a Lagrangian defined on the tangent bundle of E. 



As is well known, the solution satisfies the Euler-Lagrange equations 

dt 5p 6p ^ ^ 

In order for these equations to make sense rigorously, E must be endowed with a smooth manifold structure. This 
is achieved in [ , ]. The key concern with Problem 1.1, however, is that numerical implementation of (1.1) is 
rendered unfeasible by a lack of coordinates for the shape space E. Sadly, the solution derived via the approach 
of direct application of the variational principle lacks practical utility, despite the analytical understanding of the 
problem. The task for this paper is to investigate alternative formulations of Problem 1.1 that takes advantage of 
the analytical results in such a way that their solutions are tractable to practical implementation. 
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Introducing un-reduction Recent progress has been made for the matching problem with simple, planar, closed 
curves in the geodesic case, with i^{p, p) = Gf {p, p) where is a Riemannian metric on the shape space E = 
Emb (5^ , IR^) / Diff (5^) . This approach involved formulating a related geodesic problem on Emb (S"^ , M^) rather 
than S, [ 1 , 2]. These papers study the geodesic problem by adapting techniques commonly used for reduction by 
symmetry of variational principles, as in [ ], in a novel way. In general, given a principal G-bundle tt : Q S and 
a G-invariant Lagrangian L : TQ — > M, one may write down a reduced variational principle on TQ/G, as in [/ ]. 
This procedure is called Lagrangian reduction by symmetry. The novel insight used in [ , ] is to use this procedure in 
reverse. That is, one is interested in geodesies on the base space S, but instead of applying a variational principle 
directly, one formulates a geodesic problem on Q that reduces to the desired geodesic problem on S. That is, a 
form of un-reduction is employed, and the resulting equations on Q are more convenient to deal with numerically. 
This approach has also been applied recently in [ ] to achieve promising results for higher dimensional geodesic 
problems. 

The un-reduction procedure is distinct from the inverse process to Lagrangian reduction by symmetry, called 
reconstruction. Lagrangian reduction by symmetry relates the Euler-Lagrange (EL) equations on TQ with a new 
set of equations on TQ/G, known as the Lagrange-Poincare (LP) equations. Correspondingly, the reconstruction 
procedure relates solutions of the LP equations on TQ/G with solutions of a set of EL equations on TQ. The goal 
of un-reduction, however, is to find a parametric family of equations, called un-reduction equations, on TQ whose 
solution projects onto those of a set of EL equations onT (Q/G). Therefore, un-reduction and reconstruction only 
coincide when the component of the LP equations on T [Q/G) coincides with the desired set of EL equations. In 
fact, even when this coincidence occurs, reconstruction turns out to be a sub-manifold of an infinite dimensional 
space of un-reduction equations. Nevertheless, un-reduction equations that coincide with reconstruction do arise 
in the case of geodesic equations with horizontality conditions enforced, and variations around that theme. This 
fact has enabled the development of horizontal shooting methods for image dynamics in [1, 2]. The term 'un- 
reduction', which was coined in [ ], reflects the fact that the technique lifts EL equations from T {Q/G) to TQ, this 
mirrors the way reduction projects EL equations from TQ to TQ/G. 




Figure 1.1: A trajectory in Q from to qi that 
projects via tt to a trajectory in S from po to pi. 
Un-reduction poses a boundary value problem 
in Q such that the projected trajectories are so- 
lutions to the Euler-Lagrange equations in S. 



Some numerical difficulties persist, in that, although the geodesic in shape space may be successfully represented, 
its parametrization may still evolve in an undesirable way. For example, when the geodesic equations are imple- 
mented on Emb (5^, M^), phenomena such as clustering of data points and bad parametrizations present them- 
selves, as shown in Figure 1 of [ ]. One approach developed in [ ] and [8] deals with these issues by adding a 
step in the numerical procedure that reparametrizes the initial conditions to evolve onto a prescribed parametriza- 
tion of the target shape. This initial reparametrization procedure commutes with evolution along the matching 
trajectory, and may also iterated a number of times throughout the evolution without significantly changing the 
problem. This iterative approach effectively breaks the matching problem up into a number of sub-problems. 
While this method works well in practice, it is extrinsic in nature, and sidesteps the problem with the matching 
equations rather than dealing with it directly. It would be preferable to seek matching equations for whom the 
parametrization problem never arose. That is, for the parametrization issues to be dealt with intrinsically . Indeed, 
an intrinsic method is also described in (loc. cit.), which involves relaxing constraints on the initial conditions. The 
extra initial degrees of freedom so gained allow one to match any boundary data irrespective of parameterisation, 
however nothing is said about the parametrization dynamics along the trajectory. One of our aims in developing 
un-reduction is to provide the necessary tools to derive matching equations that intrinsically control parametriza- 
tion d5mamics along the entire trajectory. Such a derivation is given in §6, which represents, as far as we are aware. 
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the first geodesic matching equations with dynamic, intrinsic parametrization control. 

The papers reporting progress so far, such as [1, 2, 3], have constrained themselves to geodesic problems. More 
precisely, geodesies on S are obtained by projection of horizontal geodesies on Q, relative to a G-invariant Riema- 
nian metric on Q that projects to the given metric on E. Whilst the geodesic problem is important, it may be useful 
to match curves in a variety of other ways depending on the application. We shall find that in order to achieve 
the desired goals in parametrization dynamics, the geodesic properties must be partially sacrificed. Further, in 
certain situations it may be useful to assert that the curve dynamics respect different properties such as that of 
being described as a graph in the ambient space. For future work we are thinking, for example, of the geometric 
splines approach for image analysis as recently discussed in [ , ^ '^], which also fits naturally into the un-reduction 
framework that we develop here. However, applications to image analysis using geometric splines is beyond our 
present scope and will be discussed elsewhere. The framework used in developing geodesic matching needs to be 
extended in order to incorporate non-geodesic properties. Our motivation in this paper is to provide one possible 
extension by developing extra geometric tools, and to demonstrate their use on an established problem from curve 
matching. 

The purpose of this paper The objective of this paper is to investigate the extent to which the un-reduction 
approach may be used in the design of numerical methods, by investigating its geometrical context. Our investi- 
gation, inspired by the novel approach taken in [ ] that introduced reduction-by-symmetry techniques in reverse 
for the particular case of geodesic simple, closed, planar curve matching, reveals a rich geometric framework for 
un-reduction, which turns out to be broader in scope than the pioneering curve matching examples. The essential 
features of an un-reduction procedure are highlighted in purely geometric terms, giving clarity and rigorously cap- 
turing the generic notion of 'doing reduction-by symmetry in reverse'. When cast in the language of geometry, the 
un-reduction technique is sufficiently general that one may expect it to find many other productive uses, even in 
different fields such as data assimilation. See also [ ] for additional outlook toward further potential applications 
of un-reduction. 

The general theory re-applied to the curve matching problem in §6 goes some way towards answering the call for 
greater control over parametrization d5rnamics in [ ]. For example, our investigation of the geometry shows that 
the problem of quality control in parametrizations, and the extension to include potential forces may be addressed 
simultaneously . Further, the general un-reduction algorithms admit the introduction of a family of exogenous 
design factors that may be used either to formulate optimal control problems based on parametrization d5mamics, 
or for modeling purposes. This extra modeling capability may then lead to greater functionality of the solution, 
such as adaptive matching algorithms in which data points may evolve to preserve properties such as uniform 
parametrization over time. 

After introducing the required geometrical background in §2 and explaining the basic ideas in the method of 
reduction by symmetry in §3, we look at the approaches for applying the symmetry-reduction method in reverse 
as un-reduction in §4. Having achieved a good understanding of the fundamental issues, a general procedure 
for applying the un-reduction method is then developed in §5. Parallels are drawn both with particular cases of 
un-reduction and with reduction by symmetry. Finally, in §6 a new un-reduction procedure is applied to the curve 
matching problem, and comparisons are made with previous treatments. 

2 Review of geometric constructions 

This section reviews the necessary geometric tools for the formulation of un-reduction. For a more in depth 
overview and application of these tools to Lagrangian reduction, see [ ]. 

Consider a free and proper right action 

QxG^Q, {q,g)^q-9 

of a Lie group G on a manifold Q, and denote hy tt : Q Y, Q/G the associated principal G-bundle. Let g 
be the Lie algebra of the Lie group G. A connection form ^ on tt is a g-valued one-form that satisfies the following 
properties: 

1. A{q • C) = for all ^ e g and q e Q, 

2. A (vg ■ g) = Adg-iA{vq), for all e T^Q and 5 G G, 
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where 

q ■ exp(tO (2.1) 

t=0 

is the infinitesimal generator associated to the Lie algebra element ^ G g and Vq ■ g denotes the tangent lifted action 
of G on the tangent bundle TQ. 

A connection form splits the tangent bundle into vertical and horizontal distributions. That is, TQ = VQ ® HQ 
where VQ and HQ are the distributions in TQ defined fiberwise by 

VqQ = ker T^tt, H^Q = ker Aq. 

Since HqQ is complementary to keiTqir, it constitutes an Ehresmann connection on Q, and property 2 ensures the 
connection HQ is G-equivariant in the sense that Hq.gQ = HqQ-g. Such a connection is called a principal connection. 
Connection forms are in one-one relationship with principal connections. We will denote by 

-.TQ^ VQ and P'' : TQ ^ HQ 

the projections associated to the decomposition TQ — VQ ® HQ. 

The horizontal lift operator, Hor, : TpS — > HqQ, p = n (q), is by definition the inverse of the isomorphism 
Tg7r|^ Q ■ HqQ ^ TpT;. The horizontal lift operator is G-equivariant in the sense that Hor^.g [vp) — (Hor, [vp)) ■ g. 
Integrability of HQ is measured by the curvature form of A, which is defined by 

B{Vq,Uq)^ dA {Vq, Uq) + [^1 {Vq) , A {Uq)] . 

The principal connection HQ is integrable if and only if S 0. For more discussion of these points, see [. ] and 
references therein. Note that S is a horizontal form, that is B {uq, Vq) = if either Uq e VqQ or Vq e VqQ. 



2.1 Associated bundles 

An associated vector bundle is a quotient manifold of the form Q xq V := {Q x V) /G where V is a vector space 
upon which G acts linearly, and G acts on Q x F by diagonal action. The equivalence class of {q, v) G Q x V is 
denoted [g, A vector bundle structure ry : Q XqV ^ T,is given by 

TV {lq,vjv) = 7r(g) . 
The vector space operations are defined fibrewise by 

|g, vjy + A Iq, ujy ^lq,v + Xu]y , for all A e M. 

Note that q must be the same for each element representative in this relation. 

A connection form AonQ introduces a covariant derivative of curves inQ XqV which reads 

A Iq it) , V it)]y = 1(7(0, m - vit) ■ A iq{t),qmv (2-2) 
where • denotes the infinitesimal action of g on V. 

The adjoint bundle is the associated vector bundle with V = g under the adjoint action by the inverse, ^ i-> Adg-i^, 
and is denoted Ad Q := Q x a Q- The adjoint bundle is special, since it is a Lie algebra bundle. That is, each fibre is 
a Lie algebra with the Lie bracket defined by 

Sections of the adjoint bundle correspond to G-invariant vertical vector fields on Q: Consider a section ^ (p) = 
lq,£.{q)]g where tt (g) — p and ^ : Q — g is a smooth, G-equivariant map. The section ^ generates a vertical, 
G-invariant vector field € X{Q) according to 



^Q{q) ■■^q-{aq))^aq)Q{q), 
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where • denotes the infinitesimal action of g on Q, see (2.1). The properties of verticality and G-invariance are 
easily verified. 

The coadjoint bundle is the associated vector bundle with V ~ g* under the coadjoint action, /i i— > Ad*/^, and is 
denoted Ad* Q — Q x g g* . There is a natural pairing between the adjoint and coadjoint bundles given by 



By using the general formula (2.2), one obtains the following expression for the covariant derivatives of curves in 
AdQand AdQ*: 

Dtlq{t) ,fi{t)]^, ^ g(t) ,/i(t) -ad^(^(t)^^(j))pi(i) . 
Therefore, for any curve fl {t) in Ad* Q and ^ (t) in Ad Q covering the same curve p{t) e E, we have 

The connection form ^ on Q induces a map A : TQ — > Ad Q defined by 

Note, in particular, that 
for all sections ^ G T (Ad Q). 

Since the curvature form B is horizontal, it induces an Ad Q-valued two-form ;B on S defined by 

Sp{Up,Vp) = [q,Bg{Ug,Vg)]^ , Up,Vp G TpS, 

where Ug, Vq are arbitrary vectors in T^Q such that Tg{uq) = Up and Tqn{vq) — Vp. The two-form B is called the 
reduced curvature form. 

Finally, recall that the cotangent lift momentum map, J : T*Q — > g*, is defined by 

(J {aq) , Og.xB = ("9:9- Ot-qxtq 

for all ^ e g and Uq G T*Q. The cotangent lift momentum map induces a map 3 ■.T*Q Ad* Q defined by 

J(a,) = Ig,J(a,)]^.. 

3 Lagrangian reduction with a force field 

Un-reduction is closely related to classical Lagrangian reduction by symmetry. This section briefly reviews the 
main results of Lagrangian reduction by symmetry. The description given here by no means does the topic justice 
and the interested reader is referred to [7, 11, 12, 13, 14] to name just a few works on the subject. We slightly 
generalize the Lagrange reduction theorem by including the effect of an equivariant force field. Such a field is 
encoded by a fiber-preserving map F : TQ T*Q over the identity, see §7.8 in [ ]. The G-equivariance property 
reads 

F{vq ■ g) = F{vq) ■ g, for all g e G, 
where • means the tangent and cotangent lifted actions, respectively. 

Let L : TQ M be a G-invariant Lagrangian under the tangent lifted action of G on TQ and consider the associated 
reduced Lagrangian £ : TQ/G M. Fix a connection ^ on Q and consider the vector bimdle isomorphism 

a A -.TQ/G ^ TT,® AdQ over S given by 

aA{[vg]a) =T7:{vq)®A{vq). (3.1) 
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The inverse is 

{up®£) ^ [Hor^Mp + Iq (g)] ^ . 
Thus the reduced Lagrangian £ may be regarded as a map £ : TS © Ad Q M. 

Consider a G- equivariant force field F : TQ T*Q. We define the reduced force fields and F^'^ by 

F^^ ■.TJ:®AdQ^ Ad*Q, F^{p,p,d):^3{F{q,q)) (3.2) 
: TS® Ad Q^T*E, p, a), «p) q), Hor,(i;p)) , (3.3) 

where (q, g) £ TQ are such that ([g, gjc) = (p, P, ct). Note that by G-equivariance of F, the right hand side does 
not depend on the choice of (g, g) in the equivalence class. The force field F is related to F^ and i^^^as follows: 



(F(g, g), w,) = {F{q, g), Hor,(T7rK))) + ^^^(g, g), iAiw,))^ (g) 
= {Ror;{F{q,q)),Tn{w,)) + (J(F(g, g)), ^K))^.^^ 
= p, a),Tn{w,)) + (j(F(g, g)), ^K)) 

- p, a), TttK)) + {F^^ip, p, a), A{wg)) . 

It will be convenient to define also the vertical and horizontal part of F by considering the decomposition T*Q = 

{VQ)° © (HQ)", where 

VQ° T*Q I (ttg, Vq) = 0, for all Vg G K,Q} and (3.4) 

HQ° : = T*Q I (ag, Vq) = 0, for all Vq e i7,Q}, (3.5) 

are the annihilators of the vertical and horizontal distributions, respectively. We thus define 

F'' : TQ VQ° and F'" -.TQ^ HQ° 
as F'^ivq) := P^(F(w,)) e VQ° and ^^(w,) := P*{Fivg)) G and obtain the relations 

Ff^ivq)^!:* (F^'iaAiv,))) and = ^* (pA'^Ca^Ci;,))) . (3.6) 

We now state the Lagrangian reduction theorem ([11, 7, 14]) in the case that external forces are allowed. 

Theorem 3.1 (Classical Lagrangian reduction with forces) Consider a curve q{t) £ Q and define the two curves 

7r(g(t)) e S and a (t) := A{q{t)) e Ad Q, 
where A is a fixed connection omr : Q ^ Y.. Then, the following statements are equivalent: 

1. Hamilton's variational principle, 

5 r L {q{t)A{t)) dt+ r F{q{t),q{t))6q{t) dt = 0, (3.7) 
Jti Jti 

holds for variations Sq{t) ofq{t) vanishing at the endpoints. 

2. The curve q{t) satisfies the Euler-Lagrange equations with external forces: 

where Vt fs the time covariant derivative associated to an arbitrary torsion free affine connection on Q and ^ is the 
partial derivative of L relative to the same connection. 
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3. The constrained variational principle (of Lagrange-d' Alemhert type) 

£{p{t),p{t),a{t)) dt 

{{F'^'"^{p{t),p{t),a{t)),n{t)) + p{t),^{t)),5p{t))) dt = 



holds, using variations Sp{t) of p{t) vanishing at the endpoints, and variations of ait) of the form 

6<j{t) = Dtvit) - [a{t),f^{t)] - B (pit), 6p{t)) 

vanishing at the endpoints. 
4. The Lagrange-Poincare equations for £ with external forcing F^'^'^ and hold: 

(3 9) 

where ipB denotes the Ad Q-valued one-form on E defined by 

ipB:^B{p,-), 

Vt denotes the time covariant derivative relative to a torsion-free affine connection V on S, and denotes the partial 
derivative of I relative to the connection V. 



Remark. We have chosen to write the Euler-Lagrange equations (3.8) with the help of a torsion free connection V 
on Q, in order to have a global (i.e. coordinate independent) formulation. One can write as usual these equations 
locally as 

d 5L SL 

where ^ denotes the partial derivative of L in a local chart. 

Energy and momentum map. In the presence of a force field, the time derivative of the energy E{q, q) = 
g, §1^ — L{q, q) along a solution of (3.8) is 

^^Eiqit),q{t))^{Fiqit),qit)),m)- 

A direct computation, using the same arguments as in §2.7 [15], shows that in the presence of forces, Noether's 
theorem is replaced by the relation 

V (^-^) = J 9)) = J (F'^iq, 9)) ■ (3-10) 



dt \ Sq 

Euler-Poincare reduction. In the particular case Q ^ G, the G-rnvariant force field is completely determined 
by a smooth map / : fl ^ 0* and we have {F{vg),Wg) ~ {f{vgg^^), Wgg^^). In this case, the Lagrange-Porncare 
equations (3.9) recover the Euler-Poincare equations with force 

In this particular case, one observes that these equations are equivalent to Noether's formulation (3.10), since they 
can both be written as 

In the absence of forces, this recovers the usual fact that Euler-Poincare equations are a direct consequence of 
Noether's theorem. 



4 Distortions of projected dynamics 
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In order to formulate an un-reduction procedure, we must relate the Lagrange-Poincare equations (3.9) with the 
Euler-Lagrange equations (1.1) on S using the information obtained in Theorem 3.1. Here, differences between 
the Lagrange-Poincare equations and the Euler-Lagrange equations are called distortions. This terminology is 
suitable since, in the context of un-reduction, such differences are unwanted and the main goal of the un-reduction 
procedure is to eliminate them. 

The conception of certain terms in the reduced equations as unwanted barriers to the objective stands in contrast 
with applications of reduction by symmetry, where such differences instead attain physical meaning. Examples 
of this occur, for example, in [ ] in the context of molecular strand dynamics where micro-structure of a strand 
couples with the reduced filament dynamics, or again in Kaluza-Klein constructions where curvature terms are 
understood as an electromagnetic field. 

There are two distortions to contend with: 

1. Coupling distortion: One often specifies a G-invariant Lagrangian L : TQ — > M as opposed to a Lagrangian 
on TE, since coordinates on S are often unavailable. Since (3.1) shows that TQ/G = TE AdQ, the La- 
grangian L certainly depends on AdQ if it is non-degenerate. We must therefore remove any coupling 
between the Ad Q dependencies and the TE dependencies from the unreduced Euler-Lagrange equations on 
Q. In general, there is no preferred way to deal with these dependencies, and therefore this issue must be 
bourne in mind when setting up each particular un-reduction problem. This ambiguity of the Lagrangian 
on TQ is referred to hereafter as coupling distortion. 

2. Curvature distortion: The Lagrange-Poincare equations (3.9) differ from the Euler-Lagrange equations (1.1) 
in that the right hand side contains a driving term that arises from the curvature of the principal connection. 
This driving term is referred to as curvature distortion. We shall not discuss the nature of curvature distortion 
further, and point it out only for later reference when we develop methods for un-reduction in §5. 

We now describe a particular class of G-invariant Lagrangians L : TQ — > M well appropriate for the formulation 
of un-reduction. We say that the G-invariant Lagrangian decouples relative to a connection A if it can be written as 
the sum of two G-invariant Lagrangians, L — L''' + L"" , where : TQ M and L"" : TQ — > M are such that 

L'*(i;,) = L''(P'*(i;,)) and {v,) ^ {P%v,)), for all i;, e TQ, 

that is, L^{vq) depends only on the horizontal part of Vq and L'"{vq) depends only on the vertical part of Vq. Now 
we observe that the bundle isomorphism a a ■ TQ/G — > TE Ad Q induces two bundle isomorphisms 

a^lvQ/G^^Q/G^AdQ and a^l^Q/^ : i/Q/G ^ TE. 

Therefore, L'^ is completely determined by a Lagrangian £^ : TE — > M and is completely determined by a 
function : Ad Q — > M through the relations 

inorqivp)) = l^'ivp) and (^Qiq)) = i^" ([g, ^1,) . 

Note also that 

^^^*^eV^go and ^=X^eiJQ°, (4.1) 

oq op dq oa 



where the last equality implies the relation = J 



Example 4.1 (Geodesic problems) Let ^7 be a G-invariant Riemannian metric on Q, and consider the Lagrangian 
T{vq) := -^WvqWg associated to the metric. We will show that L decouples relative to a connection and we will 
compute the induced functions £^ and i^"^. 



Recall that a G-invariant Riemannian metric naturally induces a connection defined hy HqQ = ( V^Q)^ and called 
the mechanical connection. The associated connection form .4„icch is determined by the the relation 

Gq {vq,q- a) = {I{q)Anicch (vq) , <^)g.y,g, for all (J e Q, 
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where 1(g) : g — > g* is the locked inertia tensor defined by 

Qq{q-v,q- cr) = {l{q)v, fJ.)^,^,^, for all v,a eg. 

The Lagrangian L{vq) — ^ \\vq\\g obviously decouples relative to the mechanical connection, since we can write 

L K) = I \\v,\\l = \\\Pv + 1\\Ph {vq)\\l = L"{vq) + L\vq). 

In order to write explicitly, we need to introduce the vector bundle isomorphism I : Ad Q Ad* Q over the 
identity defined by 

The induced functions : Ad Q ^ M and F ■.TY.^M. are thus given by 

i^" (a) = (a) , a) and (p, p) = ^Mlg, 

where 7r*5 is the Riemannian metric on E induced hy Q. ♦ 

Distorsion in the Lagrange-Poincare equations. Example 4.1 shows that the class of geodesic Lagrangians is 
contained within the restricted class of Lagrangians on Q that decouples relative to an appropriate choice of con- 
nection form. In this case, the reduced Lagrangian e : TQ/G takes the form £ = £^ +£^<i. Consequently, 
using the relations 

S£ _ S£^ d£^'^ 6£ _ S£^ d£ _ 6£^'^ 
Sp Sp 6p ^ Sp Sp ' Sa Sa ' 

the Lagrange-Poincare equations (3.9) read: 

(4 2) 

Here the third term on the right hand side of the second equation in (4.2) is the curvature distortion, while the 
second term gives an explicit form of the coupling distortion. Note that the left hand side is the Euler-Lagrange 
operator on S for the Lagrangian £'^ : TY. — > M. Thus, the goal of the un-reduction method is to selecting solutions 
to a Lagrangian system with forcing on Q such that the right hand side of the second equation vanishes. Such 
solutions in Q project onto solutions of the Euler-Lagrange equations on E with Lagrangian £^ . 

5 Un-reduction 

Having investigated the differences between the Lagrange-Poincare equations and the Euler-Lagrange equations 
on TYj, we proceed, in light of (4.2) by using the external force to cancel the right hand side of the second 
equation in (4.2). Clearly we need to define the reduced force : TE © Ad Q — > T*E by 



_ \ ^^Ad 



F^(p,p,a):=(— ,ip^)-— . (5.1) 



For the moment we do not specify the other component : TE Ad Q Ad* Q of the reduced force. However, 
once F^'^ is fixed, then the force field F : TQ T*Q is completely determined by the equality 

Fivq) = F'^ivq) + F^vq) = vr* (F^(a^(i;,))) + A* (F^^iaAivq))) 



as was shown in (3.6). 
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5.1 The un-reduction theorem 

In order to compute the unreduced Euler-Lagrange equations associated to this forcing term F, we apply the 
variational principle (3.7). Using the same notations as in Theorem 3.1 we have 

Q = S I L''{q,q)dt + d f L''{q,q)dt+ f F'' {q, q)6qdt + F"" {q, q)dqdt 

dt 



+ / {F^{p,f>,a),5p)dt+ / {F'"'{p,p,a),A{5q))dt 
Jq Jq 

L'^{q,q)dt + {^±-^5a^+ (^^,B{p,5p) 



+ / {F^''ip,p,a),A{Sq))dt, 
Jo 

where we used (5.1) in the last equality. Now the expression for the constrained variation Sa yields 

°(.(f),*,).<.-/(|.(f)..,„).. 

since 77(i) = lq{t),A{Sq{t))l^. Therefore, the variational principle yields the Euler-Lagrange equations 

These equations may be split into horizontal and vertical parts by observing that the terms on the right hand side 
belong to HQ° and the terms on the left hand side belong to VQ° since we have 

SL'^ Jf" , Ji"" ,,,, 

= TT — ^ and — — = TT ——. (5.3) 
oq dp oq op 

Thus, the un-reduction equations (5.2) can be equivalently written as 

d 5L^ SL'' 



dt Sq 5q 



(5.4) 



where : TQ HQ° is completely determined by F^'^. Note that by the first equality in (5.3), we have the 
relation 

Note also that equation (5.5) can be rewritten as 

since F'"{q, q) € HQ°. This equation can also be deduced from the relation (3.10). 
The results obtained so far are summarized in the following theorem: 

Theorem 5.1 (Un-reduction) Let : TE R be an arbitrary Lagrangian defined on the base of a principal bundle 
TT : Q Q/G = T,. Let L = L'^ + L'" be a G-invariant Lagrangian such that L'^ o Pi^ = L'^ and L" o _P„ = L", where L'^ 
is uniquely determined by i^, and choose an arbitrary force F"" : TQ — )■ HQ°. 

Then, solutions q{t) of the equations (5.4) - (5.5) project to solutions p{t) := 7r(q(<)) of the Euler-Lagrange equations for 
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Remark. We recall here that equations (5.4) - (5.5) are the unreduced Euler-Lagrange equations on TQ for a 
Lagrangian L = L" -\- and a force field ¥ — + F^, in the special case when F*^ is constructed from the 
given Lagrangian £^ by the formula (5.1). Note that once a connection has been fixed, the choice of determines 
and vice versa. However, the choice of L" and F" is left open. The main content of un-reduction is that the 
projected djmamics on E is independent of the choice of L^' and F"" . Therefore, both these functions may be chosen 
arbitrarily, a freedom that constitutes a distinct modeling step for the particular application at hand. In practice, 
vertical d5mamics, that is equation (5.5), may be chosen to give suitable numerical properties or to add additional 
functionality to the solution being developed. An example of such a modeling procedure is given in §6. 



5.2 Relationship with vanishing momentum maps and horizontal methods for geodesic 
problems 

The horizontal methods for geodesic problems, developed and applied in [ , ], may be recovered from the general 
un-reduction procedure developed here. The horizontal approach takes the Euler-Lagrange equations on Q with 
the geodesic Lagrangian and asserts, in addition, the condition 



5L 



= 0. (5.7) 



To gain a clearer grasp of the horizontal approach for geodesic problems, consider first the more general problem 
of asserting (5.7) for a general Lagrangian L — L'^ + L" . Since 5L^ /5q G VQ°, the relation J (5L^ /6q) = is 
always satisfied, therefore (5.7) implies J {SL^ /Sq) ~ and we have 6£^'^/Sp = J {SL"" /Sq) = 0. This reduces the 
Lagrange-Poincare equations (4.2) with zero external forces to the following single equation, 

Vt = . 

6p 6p Sp 

At the un-reduced level, since we already know that SL^ /5q e HQ° by (4.1), the condition J (SL" /Sq) — implies 
that SL'" /Sq = 0. Consequently, the Euler-Lagrange equations in the vanishing momentum case take the form 

~ SLf" VL'' VL" 

— -. = . 

Sq Sq Sq 

Therefore, the vanishing momentum map condition eliminates curvature distortion, although coupling distortion 
may persist. 

Note that in general these equations do not coincide with the un-reduction equations with the vanishing momen- 
tum map condition imposed. Indeed, from (5.4) - (5.5) we obtain that these un-reduced equations read 

V,^-^.0 (5^8, 

Sq Sq 

.(|).0, (5.9, 

where the coupling distortion term VL" /Sq has now been cancelled by the horizontal forcing. 

Specializing to the particular case when L is the geodesic Lagrangian and A is the corresponding mechanical 
connection, the vanishing momentum map condition (5.7) becomes a horizontality condition. 



S_L 



= qeHQ. 



In addition, for the geodesic Lagrangian, VL" /Sq = WL/Sq = where V is the Levi-Civita connection. Thus, in the 
case of geodesic motion with vanishing momentum, the un-reduction equations are equivalent to the horizontal 
geodesic equations themselves 

Vtg = (5.10) 
= 0. (5.11) 



12 



Whilst (5.10) - (5.11) may appear over-determined, the observation that (5.10) is an Euler-Lagrange equation with 
symmetry allows the interpretation of (5.11) as Noether's Theorem. Therefore (5.10) - (5.11) are consistent, and 
(5.11) asserts a constraint on the initial conditions. Namely, A{qo) = 0, that is the initial velocity must be horizon- 
tal. Thus, un-reduction applied to the geodesic problem together with the vanishing momentum map condition 
recovers to the horizontal shooting method for geodesic problems. 

This coincidence occurs because, in addition to the property of vanishing momentum to eliminate curvature dis- 
tortion, the particular choice of the geodesic Lagrangian does not introduce any coupling distortion. Therefore, 
the forcing (both horizontal and vertical) are equal to zero, and the un-reduction procedure drops to reconstruc- 
tion of a symmetry reduced Lagrangian system. That is, the un-reduction procedure tells us to seek solutions 
of the original Euler-Lagrange equations (i.e. the geodesic equation) that satisfy the vanishing momentum map 
condition. 

As we have seen above, however, for a general Lagrangian the vanishing momentum map condition is not admis- 
sible without modification since it does not allow us to obtain the Euler-Lagrange equations for on E. In this 
case, one needs to use the un-reduction approach developed here as opposed to simply asserting the vanishing 
momentum condition. 

6 Closed, simple plane curve matching 

This section develops equations for the problem of matching closed, simple, planar curves using the un-reduction 
techniques developed in the paper. The matching problem for such curves has been treated in, for example, 
[17, 1, 18, 2]. The formulation of the problem differs between authors, however, the methods developed achieve 
consensus in agreeing that the path followed by the matching algorithms should be geodesies. In this section, 
we present an alternative perspective. Since any observer of the curve djmamics is only sensible of the geometric 
information, or 'shape', of a curve, we propose that matching requires only that shape need follow a geodesic. 
Meanwhile, the dynamics of the whole curve, which includes information about both parametrization and shape, 
may evolve in a non-geodesic fashion. 

This broader notion of curve matching allows the introduction of an entire family of matching dynamics, one of 
whom is the original geodesic d5mamics. This new family of matching dynamics, which is described efficiently by 
the un-reduction technique outlined in this paper, incorporates enough flexibility to implicitly overcome other dif- 
ficulties. Here, for example, we use the extra flexibility brought by un-reduction matching to address the unwieldy 
parametrization d5rnamics displayed by geodesic matching, as in Figure 1 in [ ], by giving dynamics that remain 
uniformly parametrized throughout the matching procedure. This achieves the objective set out in §1 to derive a 
set of matching djmamics that addressed the parametrization problems faced by curve matching algorithms in an 
intrinsic, d5rnamic way along the entire trajectory. 

Additionally, we present a member of the un-reduction family whose dynamics respect the property of curves be- 
ing written as a graph in polar coordinates, wherever it is sensible to do so. In higher dimensional problems, such 
as that studied in [ ], the graph preservation property may significantly simplify numerical implementation of the 
d5mamics. These two examples of non-geodesic matching dynamics will, hopefully, stimulate the description of 
matching d5mamics with other additional properties. 

6.1 Geometric setup 

Following [ ], consider Q = Emb"'" (S"^, M^), the smooth manifold of all positively oriented embeddings from 
to M^. Q may be thought of as the space of simple, closed, planar curves. An element c E Q contains information 
about a shape, namely a submanifold of M.'^ of S'^-type, and a parametrization of the shape. The space of shapes 
may therefore be identified with the quotient manifold S = Q/G, where G = DifF^(S'^) is the group of orientation 
preserving diffeomorphisms of 5^ which acts freely and properly from the right on Q by composition. We refer to 
[ ] for a detailed discussion of the manifold structure of Emb (S^jK^) and E, and to [17] for a discussion on the 
geometry of spaces of immersions and their quotients. 

In particular, in an appropriate topology, Emb^ (5*^,18^) is an open subset of G°° (5^ , ) . The quotient manifold, 
however, is not feasible for numerical implementation since there are no natural coordinates and analytical consid- 
erations can only be achieved via the use of equivalence classes. Therefore, it is natural to resort to an un-reduction 
procedure that takes advantage of the principal Diff^ (S'^)-bundle structure of Q over E. 
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Recall that an embedding c : 5^ — > is a smooth injective immersion that is a homeomorphism onto its image. 
Thus, Q may be identified with smooth, periodic maps Cf^ (M; M^), with, say, period 1, that are injective on each 
interval [a, 1 + a) for all a G M, and satisfy \cg \ {9) ^ for all e M. Here the subscript denotes differentiation. 

Similarly, an element f £ G := Diff^(S'^) may be identified with a smooth, strictly monotonic map / : M M 
such that / (6* + 1) = / (9) + 1 for all 9 eR. 

Since Q is an open subset of C°°(S'^, M^), a tangent vector Uc e TcQ is represented by a pair of maps (c, U), where 
U E Cj^(M, M^). The Lie algebra of G is given by the space g :— X [S^) of vector fields on 5^, which may also 
be represented by real valued, periodic functions G^ (M, M). The infinitesimal generator associated to u e g reads 

'^Q (c) = (c, UCg). 

The following are some geometric quantities which are convenient to use due to their behavior under the action 
of G. 

• Derivative along the curve 




• Length of the curve 

l{c)= [ \cg\d9 

Jo 

• Unit tangent vector 

r(c) = or r(c) = Dgc 

• Unit normal vector 

n(c) = Jt{c), where J = 

• Curvature 

/t(c) = {D0t{c)) ■ n(c) 

• Volume measure 

vol(c) = \cg\d9 

A cotangent vector in T*Q is represented by a pair of maps (c, P (g) n), where P (g) n e G^{R, (M^)*) (g) n^{S'^). 
The dual space to the Lie algebra g = X{S^) is identified with g* = X{S^)* = G^{R, R) ® Q^S^). This is seen by 
introducing the following pairings 

((c,P(g)f7),(c,t/)) := / P{U)n, and {^x®0J,u)=[ fi{u)u, 

where P®VLe T*Q, (c, U) g T^Q, fi^uj e q*, and u e q. 

The cotangent lift momentum map, J : T*Q ~¥ g*, that arises with this setup is calculated to be 

3{c,P(x)n) = {p-cg)(x)n. 

Fixing rt = vol(c) to be the volume measure on c, the momentum map becomes 

J(c,P®vol(c)) = |ce|(P-r(c))(g)vol(c) e g*. (6.1) 

The simplest metric to consider would be the metric. However, as pointed out in [ ], this metric has arbitrarily 
small geodesic distance between any two curves. The simplest metric that is of practical importance is therefore 
the metric weighted by curvature. That is, 

gc(U,V)^ f {l + K^c))(U-V)vo\{c). (6.2) 
Jo 
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For a detailed discussion of the properties of this and other metrics such as Sobolev metrics, see [17, 1]. The 
mechanical connection corresponding to G is given by 

U ■ t(c) 

Amecb{{c,U))^ (6.3) 

Indeed, since 

GciuQ{c),VQ{c)) = / (1 + K^(c)) TO|cep vol(c), 
Jo 

the locked inertia tensor is I(c)u = (l + k^(c)) \ce\'^u (g) vol(c). Therefore, the identity 

(I(c)^„,ech(c, U), u) - ((c, (7), (C, MCe)) 

implies formula (6.3). The projections associated with Amcch are 

P^'ic, U) = (c, [U ■ t{c))t{c)) and P''(c, J7) = (c, [/ - (C/ • r(c))r(c)) = (c, (J7 • n(c))n(c)). 

Thus, the horizontal-vertical split of a vector field along c is just its decomposition into normal and tangent com- 
ponents. 

Let L be the geodesic Lagrangian corresponding to (6.2), 

i(c,cO = \\\cA\l = \f^ (l + «'(c)) IqP vol(c). (6.4) 

The horizontal part of L is 

L'^icct)^ - {1 + K^{c))hic,ct f Yolic) (6.5) 
where the notation h{c, Ct) ■— Ct ■ n{c), so that P'^{c, Ct) — h{c, Ct)n{c). 

6.2 Derivation of un-reduction equations 

In order to derive the un-reduction equations (5.4) and (5.6) for the planar curve matching problem one must 
write down the Euler-Lagrange equations for the Lagrangian L'^ . This is best accomplished by using a variational 
principle. For the variational approach, one must first calculate the variations that will be used in the variational 
principle before calculating the equations. This calculation is achieved by the following lemma, after which un- 
reduction equations are derived. 

Lemma 6.1 The variations 6h, S vol(c) and 5k are 

Sh = {dt- sDg){6c-n) + Dehi6c-T) 
(5vol(c) = {Dg {dc ■ t) — {Sc ■ n) k) vo\{c), 
5k = {Dg + K^) {Sc ■ n) + DgK {5c ■ t) . 

where for simplicity we write h, n, t, k instead ofh{c, ct), n{c), t{c), k{c). 

Proof. It is convenient to calculate the commutator of differential operators, [5, Dg], before continuing. One has 

\cg\ \cg\ 
cg ■ 5cg „ 



\ce\' 



= -{T-Dg5c)De 

= {{n- 5c)k- Dg{T ■ 5c))Dg, (6.6) 
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where we used Dgr = ku. Employing (6.6), one may precalculate the variation of the tangent field, r, 

St = SDgc 

= Dg5c+[S,Dg]c 

= Dg ((r • Sc) T + {n ■ 5c) n) + {{n ■ 5c) n — Dg (t • 5c)) t 

= {{t ■ 5c) K. + Dg {n ■ 5c)) n. (6.7) 
where we used Dgn = —kt. Similarly, for the time derivative we have 

dfT = ((r ■ ct) K + Dg {n ■ ct)) n = ((r ■ ct) n + Dgh) n. 
Multiplying these variations by J we obtain for n 

5n = J5t = — (Dg {5c ■ n) + {5c ■ t) k)t and dtn = JdtT = — {Dgh + sk) t, 

where we used the notation s :— ct ■ t. 
Proceeding with the expression for 5h, one has 

5h — 5 {ct ■ n) — 5ct ■ n + ct ■ 5n 

— dt {5c ■ n) ~ 5c ■ dfji — {Dg {5c ■ n) + {5c ■ t) k) s 
^ {dt - sDg) {5c ■ n) + Dgh {5c ■ r) 

as required. 

The expression for 5vol{c) follows since 

5 vol(c) = — j — j — do 
\ce\ 

= {DgSc ■ t) vol(c) 

= {Dg {5c ■ t) — {5c ■ n) k) vol(c). 

Finally, upon noting that 5t cx Dgr cx n, and similarly reversing r and n, the derivation of the expression for 5k 
reads 

5k = 5 {DgT ■ n) 

= DgSr ■ n + [5, Dg]T ■ n + DgT ■ 5n 
= Dg {5t ■ n) + [5, Dg]T ■ n 

= Dg {{t ■5c)k + Dg (n • 5c)) + k ((n • 5c) k - Dg {t ■ 5c)) 
= {Dj + K^) {5c ■ n) + DgK {5c ■ t) , 

as required. ■ 

Proposition 6.2 The un-reduction equations for simple, planar, closed curves with the geodesic Lagrangian associated with 
Q read 



h = Dg{sh)~^^^^h^ + j^^^^{Dl{Kh')-2KhDlh) 
St = f{c,s,h) 
Ct = hn + ST. 

where f is an arbitrary smooth, G-invariant map Emb (S'\M2) x (M) x (M) (M). 
Proof. 
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The first step is to calculate the variational principle for Lagrangian (6.5) using the variations derived in Lemma 
6.1. The variational principle may be calculated as 



^ I j \{^ + ^^)^^ vol(c) dt 
Jo 

[ {1 + K^) h5hvol{c)dt+ [ [ - {1 + K^) S vo\{c) dt + [ [ nh^Sn vo\{c) dt. 
^0 Jo Jo 2 Jo Jo 

After substituting the variations from Lemma 6.1 and integrating by parts, the terms proportional to {6c ■ t) vanish 
as expected since the Lagrangian L'' is known to be degenerate. The remaining terms, those proportional to {5c ■ n), 
yield the equation 

ht = Dg {sh) + {{Dj + K^) (nh^) -h{dt^ sDe) (l + k^)) . (6.8) 

Now, replacing 5 with dt in the expression for 6k from Lemma 6.1 yields the following result 

{dt - sDe) K = {Dl + K^) h. 
Employing this relation on (6.8) one obtains 

Kh? 1 

ht = Dg {sh) - — + — ^ {{Dl + K^) {nh^) - 2Kh [Dl + k^) h) . 
Rearranging terms yields the desired result. 

The equation for c is simply a decomposition of ct into horizontal and vertical parts. That is, 

ct = P'\ct) + P^ct) = hn + ST. (6.9) 
Next, for the equation for s, note that (6.1) implies the relation 

Ai = J (^^^ = \cg\ (1 + K^) s (g) vol(c) = I(c)s (6.10) 

so equation (5.6) reads ^ (I(c)s) = J {F"" (c, Cf )) . Computing the time derivative on the left hand side and rear- 
ranging the terms yields the relation st = I(c)~^ (J {F^ (c, q)) — ^(I(c))s). Substituting in equation (6.9) on the 
right hand side yields an equation for s that may be expressed as 

St f{c,s,h) , 

for an arbitrary smooth, G-invariant map /. Of course, there is a relationship between / and J (F" (c, ct)), however, 

this relationship is not necessary for the present purposes, so it has been omitted. 

■ 

The un-reduction equations for curve matching derived in Proposition 6.2 represent a family of d5mamical equa- 
tions parametrized by an exogenous choice of smooth function / : x — >■ K. For any initial conditions 
(co, Co) e TQ and any choice of exogenous function /, the projected dynamics in Q/S, obtained by taking the 
image of a curve c, are guaranteed to follow geodesies under the metric induced on shape space by the curvature- 
weighted L'^ metric on Q. Consequently, the choice of parametrization of the boundary data in matching algo- 
rithms, and the choice of forcing may be freely chosen for convenience of the user, to pose an optimal control prob- 
lem, or to introduce extra modeling into the dynamics without needing to re-derive the equations. We present 
a couple of uses of this extra flexibility in the remainder of this section. It should be noted, however, that the 
examples we give only scratch the surface of possible uses of the extra functionality brought by un-reduction. 

6.3 Particular choices of vertical forcing and initial conditions 
6.3.1 Horizontal geodesies 



The horizontal geodesies case has already been described in some detail in §5.2. Recall that c is horizontal if s = 
(consistently with the vanishing momentum map condition (6.10)). This case corresponds to the choice f — 
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together with the initial condition sq = 0, and recovers the horizontal shooting method of [ , 2]. The equations for 
this case are given by 

C( = hn. (6-12) 

Note that by choosing / = 0, but allowing a free choice of the initial condition so, one recovers the method 
introduced in [2] to address parametrization concerns. These equations read 



h = De{s^h)-^^^^^^h' + j^^^{Dl{Kh^)~2KhDlh) (6.13) 



Ct = hn + Sqt. (6.14) 

Equations (6.13) - (6.14) incorporate an extra term sq compared with (6.11) - (6.12) which allow matching without 
the need to reparametrize the boundary data. Note that nothing is asserted about the parametrization away from 
the boundary data. Differences between the equations derived here and those in [2] are attributed to different 
approaches being used. Their equations live on X (M'^)*, the one-form densities on M^, with the equations on T*Q 
being given by a momentum map. Our equations and variational principle are formulated directly on TQ. Whilst 
the precise relation between the two sets of equations remains unclear, the method of adding an initial vector field 
on in order to enable matching parametrized boundary data is in agreement. Indeed, comparing equations 
before and after adding in the extra parameter both here and in [ ] shows that each set of equations are modified 
by adding precisely the same terms. 



6.3.2 Uniform parametric morphing 

The second case considers forcing / that ensures that the curve c is parametrized imiformly for all time. We call 
this uniformly parametrized morphing. That is, that the equation 

\ce\^l{c) / vol(c) (6.15) 
Jo 

holds for all time and for all 6 £ S^. Notice that the right hand side of (6.15) is independent of 9. [' shows 
that such curves, which are equally well described by the relation \ce\ = constant, form a submanifold U C 
Emb(S'^, M^). Differentiating (6.15) shows that a tangent vector (c, U) G TU if and only if 

\co\{De{U ■t)~{U ■n)K) = - I (JJ • n) k vol(c). (6.16) 

JO 

Indeed, if 5c — U, the variation of the left hand side of (6.15) is 

c Sc 

S\cg\ = — — — = cg ■ DeSc = Ce ■ i6T - [S,Dg]c) ^ \ce\{De (U ■ t) - {U ■ n) k) , 
\ce\ 

where we used that cg ■ dr — 0. The variation of the right hand side is 

S [ vol(c) = / {Dg {U ■t)-{U ■ n) k) vol(c) = - / {U ■ n) k vol(c) 
Jo Jo Jo 

Applying (6.16) to J7 = C( reveals the following relation 

\cg\{Dgs- hn) = - [ /iKvol(c). (6.17) 
Jo 

Solving the relation (6.17) for s in Proposition 6.2 would yield the required vertical term to keep the evolution of 
c uniformly parametrized. There is still a problem to deal with, namely that (6.16) does not respect the action of 
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G. This corresponds to breaking the G-invariance of / in the un-reduction equations, which is disallowed since 
G-invariance is critical to the derivation of the un-reduction equations. 

A solution to this issue may be found by noting that one only requires (6.16) on the submanifold of uniformly 
parametrized curves, U. The intersection of each fibre of Q — >^ E with U are the orbits of rigid rotation of S^. 
Indeed, (6.16) is invariant under the action of on Q given by (a, c{6)) i— >■ c{9 + a). Therefore, one may define / 
on U such that (6.16) is satisfied, and then extend to the whole of Q by enforcing G-invariance. 

Having described the way in which a suitable / may be found in principle, it now appears that there is a short cut 
that allows the explicit calculation of / to be b5rpassed. Note that on the submanifold U one has l/\co\ = 1. Thus 
one may multiply this factor on the left hand side only of (6.16) to obtain a relation that is precisely equivalent to the 
original on U, but is also G-invariant. This relation, applied ioU = ct, is 

l{Dgs-hK) = -[ hKvo\{c). (6.18) 



Solving (6.18) for s yields the same result as integrating the second un-reduction equation from Proposition 6.2 
with / constructed as described above. Rearranging (6.18) yields 
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Dgs = hn - - hn vol(c). (6.19) 
' Jo 

Equation (6.19) states that the rate of change of s as one moves along the curve is equal to the deviation of the 
quantity hn from its average around the curve. Since uniform parametrization is a non-local property of the curve, 
we expect that non-local terms appear in the equations, these terms take the form of the average of hn. 
Combining (6.19) with the un-reduction equations from Proposition 6.2 yields the following equations 
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s{0) = s(0) + / hn vol(c) - y / vol(c) / hn vol(c), (6.21) 
Jo ' Jo Jo 

Ct = hn + ST. (6.22) 

Equations (6.20) — (6.22) may be described as the lifting of the geodesic equations on E under the metric induced 
by G to the the submanifold, U C Emb (5^ , gf uniformly parametrized embeddings. As far as we are aware, 
(6.20) — (6.22) have not appeared in the literature. The main interest in these un-reduction equations is that they 
solve the major issue concerning the horizontal shooting method, (6.11) — (6.12). That is, the parametrization of 
solutions of (6.20) — (6.22) behave regularly for all time. This property solves the problems concerning clustering 
of data points and bad parametrizations implicitly along the entire trajectory, thereby solving the parametrization 
issues described in 51. 



6.3.3 Section morphing 

During a recent visit to London, M. Bauer and P. Harms asked whether un-reduction could be used to achieve 
geodesic morphing of shapes whilst respecting the property of being described as a graph. We worked with 
Bauer and Harms to show that this objective is indeed possible. The argument is presented here for the case of 
simple, closed planar curves. We would like to thank Bauer and Harms for their insight and collaboration on this 
point, and hope that the experience is helpful for their studies of higher dimensional matching problems, see some 
examples of which may be found in [ ] . 

To capture the setup geometrically, we begin with slight generalizations. Consider Q — Emb {Ai,J\f), which is a 
principal bundle over E = Q/G where G = Diff {M), as in [ ]. Suppose that M ~ M x 5 is a trivial fibre bundle 
over with fibre S and projection ttj^ given by projection on the first factor. A graph is then a section rj E r(7r^). 
We say a curve c{t) c Q respects the graph property if and only if c(t) c T{T:^f) n Q. That is, ttaa o c{t) = idM- 

Interesting choices could involve {Ai,Af) = (S*^, M^/{0} = x which could be used to match spheres em- 
bedded in under small deformations. The example {M,J\f) = (M^, {(x, y, z) e M^i^ > 0} = M2 x M+) could be 
used to compare different topographies. 
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Here, we shall specialize to the case of closed, simple planar curves. This corresponds to the particular choice 
{M,N) = ("S*^ , MV{0} = 51 X M+) . We shall identify with the unit circle in and fix 



:U = M.^/{{)} ^ M = S\ 7rAA(x) =x/|x|. 



Section morphing is not available for all curves in the plane, but rather is constrained to curves c e r(7r^) n Q that 
contain the origin. Considering curves up to translations in the plane, the choice of origin becomes arbitrary, and 
we shall therefore ignore the origin condition. 

The approach to deriving the un-reduction equations is similar to that taken for uniform parametrisation morphing 
in §6.3.2. That is, we derive a relation on the submanifold U :— T{t:j^) (1 Q <z Q. The submanifold U is not 
G-rnvariant, in fact, it is only preserved by the identity. Despite this, in contrast with §6.3.2, our relation turns 
out to be G-invariant without the need for modification. Therefore, the implied forcing may be extended to the 
G-invariant submanifold U ■ G c Q. This observation again allows us to identify our equations as un-reduction 
equations without explicitly deriving the vertical forcing. 

For the graph condition to remain invariant over time, the map ttj^ o c{t) must be time-invariant. This leads to the 
following calculation 

= dJ;^)^^(ct-cS-^]. (6.23) 



|c| 



From this relation we deduce that 



Ct ■ T = (c ■ t) * and h = ct ■ n = {c ■ n) * 



which results in the relation 



s = h^^ (6.24) 

Note that equation (6.24) is G-invariant, therefore it may be extended along the whole oiU ■ G. Differentiating 
(6.24) and substituting in equations from Proposition 6.2 would again result in a complicated expression for the 
exogenous parameter, /, that would be the graph preservation forcing term. Given that (6.24), we omit the calcula- 
tion. Collecting the un-reduction equations from Proposition 6.2 together with (6.24) yields the section morphing 
un-reduction equations 

ht De[h^^0)~'-^^^h^ + ^^ (6.25) 
Ct = j^^c. (6.26) 



Introducing polar coordinates in the plane, c{t, 9) = r(i, 6){cos 6, sin 6'), equations (6.25) - (6.26) become 

2rg\ K (1 + 3k^) 2 1 



ht - De[h^f}-^^^^^h^ + ^^^{Dl{^h')~2^hD^,h) (6.27) 



rt = h\ (^] +1. (6.28) 



Note that equations (6.27) - (6.28) describe a set of differential equations on functions taking values in M, whereas, 
the full un-reduction equations from Proposition 6.2 constitute a set of differential equations on functions taking 
values in M^. Therefore the particular choice of morphing that respects the section relation (6.23) has allowed us to 
integrate the un-reduction equations twice. In higher dimensional problems, this method of reduction which we 
call section morphing could prove invaluable for simplifying numerical implementation. Note that, in this example, 
clustering of data points is prevented by the graph preservation property. The only parametrization effect that may 
cause a problem is the distance between data points becomes large as the radius maxg r{6, t) increases. This issue 
may be averted by first scaling the shapes to be matched, and then using section morphing. A detailed discussion 
on such developments is beyond the scope of this paper, and is left to future work. 
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